Quantum dimer models and effective Hamiltonians on the pyrochlore lattice 
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We study a large- JV deformation of the 5 = 1/2 pyrochlore Heisenberg antiferromagnet which 
leads to a soluble quantum dimer model at leading non-trivial order. In this limit, the ground 
state manifold - while extensively degenerate - breaks the inversion symmetry of the lattice, which 
implies a finite temperature Ising transition without translational symmetry breaking. At lower 
temperatures and further in the 1/N expansion, we discuss an effective Hamiltonian within the 
degenerate manifold, which has a transparent physical interpretation as representing dimer potential 
energies. We find mean-field ground states of the effective Hamiltonian which exhibit translational 
symmetry breaking. The entire scenario offers a new perspective on previous treatments of the 
SU(2) problem not controlled by a small parameter, in particular showing that a mean-field state 
considered previously encodes the physics of a maximally flippable dimer configuration. We also 
comment on the difficulties of extending our results to the SU(2) case, and note implications for 
classical dimer models. 

PACS numbers: PACS numbers: 75.10.Jm, 74.20.Mn 71.10.-w 



I. INTRODUCTION 

One of the central open questions in the study of frus- 
trated systems is what happens when a classically highly 
degenerate magnet is subjected to violent quantum fluc- 
tuations. A model system which has played an impor- 
tant role in the discussion of this question is the nearest 
neighbour spin 1/2 Heisenberg antiferromagnet on the 
pyrochlore lattice. This lattice consists of corner-sharing 
tetrahedra and exhibits a massive classical degeneracyi 

Very little is known reliably about the quantum model: 
exact solutions are unavailable in d = 3 + 1 , Monte Carlo 
simulations are frustrated by the sign problem, and the 
pyrochlore lattice - being three dimensional and having 
a unit cell of four spins - does not yet lend itself to exact 
diagonalisations . 

In an important, if somewhat cryptic, initial piece of 
analytical work back in 1992, 2 Harris, Bcrlinsky and 
Bruder (HBB) took the approach of considering the 
bonds belonging to one of the two sublattices - "weak 
bonds" - of tetrahedra perturbatively (the tetrahedra of 
the pyrochlore lattice are arranged on the bipartite dia- 
mond lattice, which in the current context can usefully 
be thought of as two interpenetrating face-centred cu- 
bic lattices, see Fig. The tetrahedra on the other 
sublattice, where the "strong bonds" reside, are thus ini- 
tially decoupled, with doubly degenerate dimerised sin- 
glet ground states, parametrised by a pseudospin a. The 
basic idea of HBB is then to switch on the weak bonds 
perturbatively, thus generating an effective interaction 
between the as on different sites. This interaction deter- 
mines the eventual ground state, and in itself defines a 
difficult quantum problem. 

A number of authors have since developed ideas based 
on expansions for pseudospins on an fee lattice^*^ 



Tsunetsugu's detailed work extended the scope of the 
mean-field theory^ while Berg, Altman and Auerbach 
implemented a sophisticated numerical procedure pro- 
jecting the Hamiltonian onto this Hilbert spaced 

These approaches have in common that their starting 
point has a lower symmetry than the initial Hamiltonian, 
as a distinction between the two sublattices of tetrahe- 
dra discards the inversion symmetry of the pyrochlore 
lattice. In the same spirit, a perturbation theory in the 
weak bonds but for the full density matrix, which unlike 
the previous papers does not impose a restriction to the 
dimerised singlet subspace, at any rate found spin order- 
ing to be absent but did not study the presence of bond 
ordering^ 

In this paper, we take a different route by considering 
a Hamiltonian with an enlarged internal, rather than re- 
duced point, symmetry: we study the spin-1/2 problem 
on the pyrochlore lattice via a large- V quantum dimer 
model approach. This approach does not expressly break 
the symmetry between the two sublattices of tetrahedra. 
It provides an (artificial) small parameter, 1/N, which 
we will use to obtain an analytical solution to the first 
non-trivial order of the quantum dimer model. 

Very unusually, we find that the quantum dimer model 
generated at 0(1/ N) is exactly soluble^ There turns out 
to be a phase transition in the Ising univsersality class at 
a temperature T ~ 0(1/ N), where the inversion symme- 
try of the pyrochlore lattice is spontaneously discarded. 
The set of states selected at this order preserves an ex- 
tensive entropy. 

The question what happens within this subspace at 
higher order in 1 /N is closely analogous to that posed in 
Refs.0@0H as the respective starting Hilbert spaces are 
isomorphic. Our approach thus provides a way of reach- 
ing this starting point dynamically, i.e. by spontaneous 
rather than explicit symmetry breaking. 
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We then s how that the effective Hamiltonians obtained 
in Refs. I2ll3l l4l5l belong to a family which has a simple 
interpretation as dimer potential terms in such an ap- 
proach; as the natural basis beyond 0(1/N) is not a 
dimer one, however, they define a quantum Hamiltonian. 

These Hamiltonians have been studied using different 
approaches We show that a supertetrahedral or- 
dering pattern proposed by HBB and Tsunetsugu can be 
represented by a simple maximally flippable dimer config- 
uration. We also show that, within the framework of the 
mean-field theory, there arc further configurations close 
by in energy, leading to numerically small characteristic 
energies. This implies that if there is translational sym- 
metry breaking, the final ordering pattern might not be 
assumed until a temperature low compared to finite-size 
gaps in exact diagonalisations or to other perturbations 
in real compounds. 

Our results thus provide an intuitive picture of the 
dominant physics captured in a class of theories pro- 
posed for the pyrochlore quantum antiferromagnet. As a 
byproduct, our results also imply that a classical dimer 
model on the pyrochlore lattice with a simple potential 
term will have an inversion and translation symmetry 
broken ground state with residual entropy. 

In the remainder of this paper, we first introduce the 
Sp(iV) quantum dimer model (Sect. ITTf> , followed by its 
exact solution (Sect. IIII(I . In Sect. IIVI we discuss the 
structure of the problem at higher orders in 1/N and 
in Sect. Ivl make a connection to previous work via the 
idea of maximally flippable dimer configurations. After 
some remarks on further extensions (higher-order loops, 
other lattices and SU(2) spins, Sects. IWVTTTIl . we close 
with a short discussion of further questions raised by this 
approach. 



II. THE Sp(iV) DIMER MODEL 

An Sp(Af) dimer model is obtained from an SU(2) 
spin problem via enlarging the symmetry of spin space 
from SU(2)^Sp(l) to Sp(iV). Rewriting the Heisenberg 
Hamiltonian in terms of operators sj-, which create a 
singlet on the bond between spins i and j, gives 

m K } 

where N = 1 for SU(2)~Sp(l). The same form of H 
holds for the generalisation to Sp(7V), with s\j now cre- 
ating Sp(A^) singlets. Details for this and the following 
steps are provided in the Appendix. 

A. The dimer model at leading order 

A hardcore dimer model is obtained from this Hamilto- 
nian at leading order— as one takes N — > oo— At this or- 



der, any nearest-neighbour dimer covering is an orthogo- 
nal, degenerate eigenfunction of the Sp(AT) Hamiltonian, 
with the ground state energy per site of Eq = — J/2. 

The ground state is thus exponentially degenerate in 
the volume of the system, although the precise value 
of the degeneracy for three dimensional lattices is not 
known. As the temperature is lowered below T ~ 1, 
the Sp(A^) magnet enters the dimer manifold of ground 
states. It is possible in principle that this manifold al- 
ready incorporates some form of order, in which case this 
would happen via a phase transition. However, it appears 
more likely that the correlations averaged over the dimer 
manifold remain short-ranged for the case of the non- 
bipartite pyrochlore lattice^ in which case the restric- 
tion to the dimer manifold has the form of a crossover. 

The dimer wavefunctions obtained here are isotropic 
in spin space. This means that any further symmetry 
breaking occurring at subleading order can only be of a 
spatial variety. 

B. Derivation of the quantum dimer model 

The degeneracy between different dimer coverings ex- 
ists only at leading order; at higher order, the model ac- 
quires a non-trivial quantum dynamics. This dynamics 
is determined along the lines pioneered by Rokhsar and 
Kivelson in the context of the SU(2) Heisenberg models 
They derived a Hamiltonian in the space of dimer cover- 
ings by formally carrying out an expansion in a parame- 
ter x which in fact has the finite value of 1/2, as we will 
briefly describe belowiii*i£ In the Sp(A r ) model, this pa- 
rameter is (artificially but) truly small: e = 1/2N. The 
derivation of the dimer Hamiltonian is hence completely 
analogous, albeit rigorously organised, in the orders of 
the small parameter. Details are again given in the ap- 
pendix. 

At 0{e 1 ), the resulting quantum dimer Hamiltonian is 
very simple. It reads: 

Hqdm = 2Jen . (2.2) 

Here, □ stands for a resonance term around a closed loop 
of length four, i.e. a kinetic term which exchanges oc- 
cupied and empty links if they alternate around such a 
loop. Note that this quantum dimer model, unlike the 
Rokhsar-Kivelson one^i only contains a kinetic term— 

III. EXACT SOLUTION OF THE QDM: 
PARTIAL ORDER BY DISORDER 

The quantum dynamics induced by Eq. 12.21 can only 
lead to resonances between two dimer configurations 
which can be transformed into one another by moving 
exactly two dimers. Such moves are only possible for two 
dimers on a single tetrahedron, as the shortest closed 
loop not confined to a single tetrahedron has length six 
and therefore would require moving three dimers in order 
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FIG. 1: The pyrochlore lattice, a network of corner-sharing 
tetrahedra. The tetrahedra on the two ('up' and 'down') sub- 
lattices are shaded differently. The centres of the tetrahedra 
of either sublattice form a face-centred cubic lattice. 

to satisfy the hardcore constraint in both the initial and 
final configurations. On a single tetrahedron, there are 
three possible dimerisations in Sp(iV), and the Hamilto- 
nian matrix elements between them read, at this order: 

/0 1 IX 
Htet = 2Je 1 1 . (3.1) 

Vi i o/ 

This matrix has one non-degenerate eigenvalue 4 Je with 
corresponding eigenvector | 4>a) = (1, 1, l)/v / 3- The other 
two eigenvalues are degenerate at — 2Je and eigenvec- 
tors = (l,x,x 2 )/V3 and \<j>J) = (l, X 2 ,x)/V3, with 
X = exp(27ri/3). These two eigenvectors form an E rep- 
resentation of the tetrahedral group Td- Note that the 
asymmetry of the spectrum under J <-> — J is a mani- 
festation of our inability - in contrast to the case of the 
squareii and triangular— lattices - to choose the sign of 
the overlap matrix elements at will. 

A single tetrahedron at this order is thus occupied by 
or 1 dimers and gains no energy from resonance moves, 
or it is occupied by 2 dimers, in which case it gains an 
energy — 2 Je < 0. We thus need to maximise the number 
of tetrahedra occupied by two dimers. 

As the number of dimers equals the number of tetra- 
hedra, this is done by putting no dimers on one half the 
tetrahedra, and two each on the remaining tetrahedra. 
As two tetrahedra with two dimers each cannot be neigh- 
bours, and as the tetrahedra reside on the bipartite di- 
amond lattice, this implies that one sublattice of tetra- 
hedra, containing the 'down' tetrahedra, say, is empty, 
whereas the 'up' tetrahedra have two dimers each. This 
gains an energy of —2Je per up tetrahedron, and hence 
—Je/2 per site. 

The quantum dimer model up to this order thus spon- 
taneously breaks the inversion symmetry through the 
sites of the pyrochlore lattice (this operation exchanges 
up and down tetrahedra); the point group is broken down 
from Oh to Td. This is, to our knowledge, the first model 
where such a symmetry breaking can be demonstrated 
to occur; this feature is not unimportant as most treat- 
ments use this symmetry breaking as part of the starting 
assumptions. We also note that this is a rare case where a 
quantum dimer model is analytically soluble - normally, 
finding the solution is a hard problem requiring intensive 
numerics in itself. 7 



Despite this symmetry breaking, the ground state at 
0(e 1 ) retains a large degeneracy as each single tetra- 
hedron continues to exhibit a two-fold degeneracy; the 
residual entropy per spin is hence S = (fcs/4) In 2, where 
ks is the Boltzmann constant. 

We thus find a transition at which the point group 
is reduced by discarding inversion symmetry, but spin 
rotational as well as real space translational symmetries 
remain intact. This is an example of (only partial) order 
by disorder in a quantum frustrated system. 



IV. THE FATE OF THE RESIDUAL 
DEGENERACY 

So far, to 0(e), we have been fortunate in dealing 
with a quantum dimer model with a simple structure. 
Configurations differing in the distribution of dimers be- 
tween tetrahedra define disconnected dynamical sectors, 
the ground-state energies of which are given by the num- 
ber of doubly-occupied tetrahedra, with the correspond- 
ing wavefunctions being outer products over tetrahedron 
wavefunctions. Note that such a simple structure is 
generically absent for other lattices, the leading dimer 
model typically not being exactly soluble. 

Here, these problems are deferred to the next or- 
der, e 2 , where the degrees of freedom are given by one 
pseudospins-1/2 for each doubly occupied tetrahedron. 
For the ground-state sector, these pseudospins reside 
on the face-centred cubic lattice defined by one sublat- 
tice of tetrahedra. These pseudospins encode the E- 
representation of Td provided by the degenerate wave- 
functions \4>±). 

Nonetheless, further insight can be gained by contin- 
uing to use a dimer basis, on the understanding that a 
diagonal operator in the dimer basis will not be diag- 
onal in the pseudospin basis. This approach will lead 
us to a simple interpretation of the effective Hamiltoni- 
ans of Refs. QEatlOjl Such a connection is possible as, 
formally, the pseudospins on a face-centred cubic lat- 
tice are also the starting point of those studies, where 
the three (linearly dependent) SU(2) dimer coverings of 
an isolated tetrahedron reduce to a two-dimensional E- 
representation. 



A. Effective Hamiltonians 

Let us consider the possible dimer operators we can 
write down beyond the ones for the loops of length 4 al- 
ready included at 0(e). The next simple kinetic term 
is that for a hexagon, O, which is in principle generated 
at 0(e 2 ). However, such a term moves dimers from the 
occupied sublattice of tetrahedra to the unoccupied one, 
and thus does not have a matrix element of 0(1) within 
the pseudospin subspace in which we are doing degen- 
erate perturbation theory. It will, however, contribute 
via a 'virtual' process by first shifting three dimers onto 
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the empty tetrahedra and then back again. This process 
is diagonal in dimer basis, pictorially represented by a 
potential term of the form 

o 

where the symbol stands for an operator which 

projects out the states with zero amplitude of dimer oc- 
cupancy on the three bonds of the occupied tetrahedra 
of a given hexagonal loop. Here, the symbol implies 
a sum over all such hexagonal loops. 

The physics of this term is most transparent in the 
dimer basis. However, it is off-diagonal in the basis of 
states we are doing perturbation theory in. Let us there- 
fore rewrite the dimer potential term in pseudospin basis. 

To do this, we write the wavefunction on each occupied 
tetrahedron in terms of two angles, a and 0: 

\9,a) = sin(a/2) exp(i0/2) \(f> + ) 

+ cos(a/2)exp(-£0/2) |<^) . (4.2) 

Next, we observe that the probability of finding a dimer 
on a given link of the tetrahedron is given by 

[1 + sinai cos(0 4 - u 0ji )] /3 , (4.3) 

where i/ 0) i — 0,±27r/3 according to which of the three 
possible dimer pairings is selected. 

Using this parametrisation, the potential term for the 
pseudospins on the face-centred cubic lattice has the form 

3 

H fcc = - n t 1 + sin a * c ° s ^ - / 3 • ( 4 - 4 ) 

O i=l 

Here, the angles v 0i i depend on the loop under consider- 
ation and are such that 9 = v Q ^ together with a = ir/2 
corresponds to the state which maximises the dimer am- 
plitude on the link which forms part of the hexagonal 
loop. Put another way, the factors in the product are 
the probabilities of finding dimers on the respective links. 
The sum on i runs over the tetrahedra to which these 
links belong. 

This form can in turn be rewritten by introducing pseu- 
dospins (represented by Pauli matrices a) on each tetra- 
hedron with the quantisation axis along the azimuthal 
axis a = 0, and e Q) j denoting the preferred directions 
9 = v 0t i and a — ir/2: 

3 

H /-=-En[( i+T ) +o! -^]/ 3 - ( 4 - 5 ) 

O i=l 

For the hexagonal potential, Hq, one has T = 0. T ^ 
corresponds to an additional potential term involving 
only two dimers on a hexagon, pictorially represented as 

nn- 

The form of Hf cc is in fact the same as that of the 
effective Hamiltonians obtained in Refs. 2.3 4§5j. The in- 
termediate steps in the algebra differ between Sp(iV) and 



SU(2) on account of the non-orthogonality of the dimeri- 
sations of a tetrahedron in SU(2) [see Sect. IVlIj . Most 
importantly, the v Q ^ are shifted by tt. 

However, in the end the effective SU(2) Hamiltoni- 
ans correspond to dimer potentials like Hq (Eq. 14.1(1 
in the same way that Hf cc does. This common equiv- 
alence states in a crisp geometrical form that both the 
SpCi VHEq.lOt and the SU(2) effective Hamiltonians of 
Refs. r2l3l l lll5l encode dimer potentials. 



B. Mean-field theory for Hj cc 

Even with the restriction to this class of Hamiltonians, 
one is still left with the hard problem of minimising Hf cc , 
which has the form of a S = 1/2 Heisenberg model with 
spin-orbit coupling. 

A proper quantum mechanical treatment of this prob- 
lem is beyond the scope of this paper. Rather, to shed 
some light on the physics unearthed in previous studies, 
we considered this problem in a mean-field theory, treat- 
ing the pseudospins- 1/2 as being effectively classicaliS*^ 

For T = 0, the best state we have discovered so far has 
a four sublattice structure, as do the states of Refs. Ul 
HI with angles a = n/2 and {(9;} = {0, 0, tt - 5, n + 5}, 
where S = arccos(7/8) « 0.5054. 

This state has an energy per loop (and therefore per 
spin), A, with A = -289/3456 « -0.08362. The highest 
possible energy, realised for 9 L = tt, gives A = 0. Even 
a state with S significantly different can have a closeby 
energy: for example, for 6 = 0, one obtains A = —1/12 w 
—0.08333. This state has a higher symmetry as it only 
represents a two-sublattice ordering. 

As T is changed, S evolves, and it reaches 5 = ir/3 
at T = —1/2; at this point, it in fact becomes possible 
to disorder one of the four sublattices at zero cost in 
energy; this we call Harris-Bruder-Berlinsky-Tsunetsugu 
(HBB-T) state. For T = 0, this configuration {0J = 
{0, 02, 27r/3, — 27r/3} is no longer optimal in our case but 
still rather close: A = -35/432 w -0.08102. 

These results raise two points. The first is the observa- 
tion that the particular four-sublattice ordering with one 
disordered sublattice is not entirely generic, as even in 
mean-field theory, it does not correspond to an extended 
parameter region in T. 

The second is the recurring theme of the emergence 
of a numerically small energy scale. There are configura- 
tions the symmetry properties of which differ significantly 
but whose energies are very similar. The configurations 
discussed above are degenerate at O(e ) and 0(e 1 ), and 
their variational energies with respect to Hf cc differ only 
by a few percent. 

These facts suggest that the magnet will break transla- 
tional symmetry - by assuming whatever ordered config- 
uration it chooses - only at very low temperatures, pro- 
vided quantum fluctuations do not prevent this ordering 
altogether. 
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V. THE HBB-T STATE AS A DIMER STATE 

Given the Hamiltonian Hf cc is most simply written in 
dimer basis, it is natural to ask whether the mean-field 
ground states have a natural interpretation in terms of 
dimer coverings. This is most pertinent for states with 
angles a = ir/2 and 9 integer multiples 2n/3, as these 
correspond to maximal dimer amplitudes on one of the 
three pairs of links^ and are thus naturally identified 
with a dimer configuration. 

The HBB-T state in fact has a very simple such inter- 
pretation: it can be represented as the dimer configura- 
tion with the maximal number of flippable hexagons, that 
is to say as the ground state of the classical dimer model 
with a hexagonal potential of the type Ho [Eq. (|4.1|) ]. 
Why one sublattice of tetrahedra remains disordered for 
the maximally flippable configuration also becomes ap- 
parent in this language. 

As an aside, we note that the appearance of such maxi- 
mally flippable configurations is a generic feature of quan- 
tum frustrated systems. In a nutshell, a large number of 
flippable loops implies a wide range of possible fluctua- 
tions, and hence undetermined degrees of freedom. These 
the perturbation can make use of by arranging them to 
its liking. Such a route to quantum order by disorder 14 
has been discussed in detail in the context of frustrated 
transverse field Ising models. 15 

To demonstrate this, we first need to show which clas- 
sical dimer states optimise this problem; then, we render 
the HBB-T state in a way which makes its identity man- 
ifest. 



A. Maximally flippable configurations 

f I 

To maximise the number of flippable hexagons, 71q, 
we note that the hexagons lie in the [111] kagome planes, 
and that each bond of the pyrochlore lattice is part of 
hexagons in two different planes, e.g. [Ill] and [-1-11]. 
Thus, each dimer can be part of at most two flippable 
hexagons; as the number of hexagons, no, equals the 
number of spins, and hence twice the number of dimers, 
and as three dimers are needed for a flippable hexagon, 
this provides an upper limit of nil/no = 1/3. 

However, it turns out to be impossible to saturate this 
bound, as can be seen by explicitly constructing con- 
figurations which are locally maximally flippable. Con- 
sider first a configuration with independently flippable 
loops. In this case, one needs at least three dimers 
per resonant loop, and one obtains for the upper bound 
Uq/uo = 1/6. A more favourable state would be one 
in which two flippable hexagons share one dimer in two 
different kagome planes, e.g. [Ill] and [-111]- For such a 
local configuration one needs at least five dimers for each 
pair of resonant loops, and this yields the upper bound 

n f o/ n O = 1/5- 

This local configuration might in principle be used as 




FIG. 2: A flippable supertetrahedron consists of six dimers, 
which form four flippable hexagons in the four {111} kagome 
planes, and four surrounding dimerfree tetrahedra of the 
down-sublattice, such as (1234). 

a building block for a new configuration if other dimers 
are added locally to form other flippable hexagons. There 
is only one possibility of adding a sixth dimer to build 
a third flippable hexagon, and this automatically yields 
a forth one, with one flippable hexagon in each kagome 
plane. This local configuration of four flippable hexagons, 
constructed by six dimers, has the form of a flippable 
supertetrahedron shown in Fig. [5] 

This flippable supertetrahedron has four dimerfree 
tetrahedra at its corners. None of the dimers emanat- 
ing from the outer vertices of these tetrahedra can be 
part of a flippable hexagon, and thus at least 2 dimers 
are lost (four dimers each shared between two flippable 
supertetrahedra) . 

To obtain four flippable hexagons, one thus needs to 
invest 6 + 2 = 8 dimers, and hence 16 sites. As the 
number of sites equals the number of hexagonal loops, 
this establishes the upper bound of 

n f o/no < 1/4 • 

A consequence of the local flippable supertetrahe- 
dron structure is that an optimal configuration, with 
n o/ n O — 1/4 automatically breaks the inversion sym- 
metry, i.e. all up-tetrahedra contain two dimers and all 
down-tetrahedra arc dimerfree. 

As an aside, we note that a classical attractive poten- 
tial for dimer pairs on a tetrahedron would have selected 
the dimer manifold (with entropy ks In 3 /4 per site) with 
all dimers on one sublattice of tetrahedra. The attractive 
hexagonal loop term (by itself or in combination with this 
term) in addition imposes the creation of the supertetra- 
hedra. 



B. Connection to the HBB-T state 

An optimal flippable supertetrahedron tiling is shown 
in Fig. |3| The dimers which do not take part in the 
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FIG. 3: Optimal dimer configuration, (a) Gray tetrahe- 
dra form the up-sublattice, and white tetrahedra (dimer- 
free) form the down-sublattice. Four flippable supertetrahe- 
dra meet in one ?-tetrahedron of the up-sublattice. The dimer 
configurations of these ?-tetrahedra may be chosen arbitrar- 
ily. The flippable supertetrahedron on the right hand side is 
represented in such a manner to make clear the relation to 
the mean-field configuration of Ref . 0. (b) Cubic lattice cell 
of the flippable supertetrahedron configuration. The flippable 
supertetrahedron is shown in red. 



flippable supertetrahedron formation are found on the 
?-tetrahedra in Fig. [3^, and may be chosen freely. The 
ground-state degeneracy of the optimal configuration 
therefore remains macroscopic, and the entropy per site is 
S e 2 = (fce/16)ln2.3 Because the up-tetrahedra form an 
fcc-lattice, the ?-tetrahedra now live on a cubic lattice 
and connect the flippable supertetrahedra. This is remi- 
niscent of the HBB-T state; indeed, the cartoon for this 
state is precisely our state of flippable supertetrahedra, 
see Fig. [3J where we have shown the flippable superte- 




FIG. 4: Twelve-bond loop connecting two nearest neighbours 
of ?-tetrahedra (blue dashed-dotted line). The dimerisation 
of the connected ?-tetrahedra has to be the same (here: blue). 



trahedra, together with the dimer configurations on the 
shaded tetrahedra as dictated by the HBB-T state. How- 
ever, by focussing on the up-tetrahedra, rather than the 
flippable loops, the geometry is harder to visualise. 

One may ask whether there exist other tilings of flip- 
pable supertetrahedrawhich are also optimal. Indeed, if 
one chooses the ?-tetrahedron in the cubic lattice cell in 
Fig. 03 in such a manner that its dimerisation coincides 
with that of the tetrahedron at lattice point [110] (3- 
dimerisation) , one may displace one layer of cubes with 
respect to a neighbouring layer in the [110] direction, as 
indicated by the gray arrow. 

However, this one choice forces us to fix the dimerisa- 
tion of all ?-tetrahedra at the interface of the displace- 
ment. The resulting interface energy will thus generically 
win over the displacement entropy and lift this degener- 
acy. 



VI. HIGHER-ORDER LOOPS 

In the absence of an exact solution at 0(e 2 ), our 
present approach does not per se justify pursuing the 
physics of longer loops. For completeness, however, we 
mention some geometric facts about the entropies of the 
residual degrees of freedom in the dimer model with the 
uniform four-sublattice ordered state, with one disor- 
dered sublattice, as a starting point. We will discuss 
this topic more completely separately^ 

The as yet undetermined dimers can only participate 
in virtual resonance loops of length 12 (involving mov- 
ing six dimers). One, rather contorted loop, can con- 
nect a pair of nearest-neighbour ?-tetrahedra (Fig. 0}. 
The other involves three ?-tetrahedra which are mutually 
next-nearest neighbours. This is displayed in Fig. 

The nearest-neighbour loop can only (virtually) res- 
onate if the dimers are oriented perpendicular to the 
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FIG. 5: Next-nearest-neighbour loops, (a) Closest packing 
of 12-bond loops connecting next-nearest neighbours of ?- 
tetrahedra (dark gray). The tetrahedron in light gray is 
part of the other square sublattice. The superstructure is 
again an flippable supertetrahedron, but with unequal bonds 
(flippable (super) 2 tetrahedron). (b) Elementary lattice cell 
of the flippable (super) 2 tetrahedron configuration. The ?- 
tetrahedra on the cubic sublattice, which are connected by 
next-nearest-neighbour loops, are represented by full circles 
(black and coloured). The open circles represent ?-tetrahedra 
on the second cubic sublattice. The four resonant next- 
nearest-neighbour loops are those connecting AFD, BED, 
BFC, and AEC. 



line joining the ?-tetrahedra, and it thus favours the for- 
mation of rods of ?-tetrahedra with identical alignment, 
somewhat reminiscent of the supertetrahedral ordering 
pattern found in Ref. 0. 

The next-nearest neighbour loops only connect tetra- 
hedra living on the same sublattice of the (bipartite) 
cubic lattice on which the ?-tetrahedra reside. Trying 
to maximise the number of flippable next-nearest neigh- 
bour loops thus leads to two identical copies of the same 



problem. Amusingly, the maximally flippable n.n.n.- 
loop configuration is self-similar: one now finds flippable 
(super) 2 tetrahedra, with one out of four ?-tetrahedra re- 
maining again undetermined (Fig.|^t>). 

The maximal total number of flippable loops of length 
12 is obtained by taking the two copies of this flippable 
(super) 2 tetrahedron state and arranging the ?-tetrahcdra 
to satisfy the rod-arrangement dictated by the nearest 
neighbour loops. 



VII. THE SU(2) LIMIT 

The Rokhsar-Kivelson quantum dimer model for SU(2) 
spins is based on an expansion in the same small param- 
eter, e, appearing in the overlap matrix. However, this 
parameter is not in fact infinitesimal - 1/2, in fact - and 
this leads to a certain arbitrariness in assigning an order 
to a given matrix element as a change in its order can 
be offset by a numerical factor 2. This is most plainly 
seen for the case of the constant 1/47V 2 ~ 0(e 2 ) in the 
Hamiltonian, which in the case of SU(2) simply becomes 
a constant 1/4 ~ O(e ). Here, we evaluate the matrix 
elements between two dimer configurations exactly and 
assign them the order of the overlap between the two 
configurations. 

On the pyrochlore lattice, the main obstacle lies in the 
fact that SU(2) dimerisations of a tetrahedron are not 
linearly independent. In the quantum dimer model, this 
shows up in the form of a non-invertible overlap matrix 
S. In such a situation, Eo . IA5I is no longer meaningful. 

One can adopt two different stances with respect to 
this observation. Firstly, one can hope that - by conti- 
nuity - the large- N physics will nonetheless be relevant 
to the SU(2)~Sp(A = 1) case. Alternatively, one can de- 
cide that SU(2) is special, and that the physics discussed 
so far in this paper will qualitatively have to be modified 
for the N = 1 case. 

Here, it is useful to note that one also finds that the res- 
onance energy of a single tetrahedron vanishes - for the 
SU(2) RK model, the matrix elements in Ea. 13. H are zero. 
Let us thus sidestep the orthogonality issue by ignoring 
the length four loops altogether, with a leading term of 
the QDM proportional to the hexagonal resonance term 
— Je 2 0. This is fundamentally different from Hf cc in 
that it represents a kinetic, and not a potential, dimer 
term. The ground state of this quantum dimer model 
(and whether it breaks the lattice inversion, or indeed 
any other, symmetry) is not known. Candidate varia- 
tional states for this problem can be constructed along 
the lines explored in Ref. yjj One such candidate is the 
maximally flippable supertetrahedron-configuration dis- 
cussed above. 

As the flippable hexagons there are not independently 
flippable, other candidate states can be constructed to 
maximise the number of maximally flippable hexagons; 
one example can be obtained from the proposal of Refill^. 
This has one in six hexagons flippable and does not break 
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the inversion symmetry. We leave a detailed analysis of 
this quantum dimer model as a subject for future study. 



VIII. OTHER LATTICES 

The approach developed here can easily be extended to 
other lattices. However, the exact solubility discovered 
here is in general not encountered elsewhere. In fact, the 
situation can differ from the pyrochlore one in several 
respects 

Firstly, the lattice under consideration may not admit 
any dimer coverings. Secondly, all dimer states may have 
a fixed number of flippable shortest loops of even length, 
thereby rendering the leading order dimer model triv- 
ial. Thirdly, the shortest loops may be inequivalent or 
overlapping in a way so as to destroy the solubility. In 
addition, the appropriate low-energy sector may at any 
rate be better described by a different type of effective 
models 



IX. CONCLUSION 

In conclusion, we have presented an exactly soluble 
deformation of the highly frustrated quantum pyrochlore 
antifcrromagnet. To our knowledge, this provides the 
first instance in which the inversion symmetry of the py- 
rochlore lattice by itself is spontaneously broken. 

From an experimental perspective, this is perhaps the 
most interesting observation in this paper: upon lower- 
ing the temperature from the paramagnetic phase, the 
current scenario implies the presence of an initial Ising 
transition into a non-magnetic phase with finite entropy 
and full translational and rotational symmetries. 

Of course, the presence of such a symmetry broken 
phase is in a se nse already implicit in the starting points 
of Refs. ElFil l llol s With respect to these, our approach has 
provided a simple physical picture by shedding some light 
on the 'natural' degrees of freedom arising in such treat- 
ments, the hexagonal resonating loops. This reinforces 
the idea of an non-magnetic ordering pattern with an 
enlarged unit cell. 

Our work therefore goes some way towards provid- 
ing a rationale for the sequences of symmetry-breakings 
discussed beforem&£& Whether or not the scenario dis- 
cussed here provides the appropriate framework for un- 
derstanding the nearest neighbour S=l/2 pyrochlore 
Heisenberg antiferromagnet remains an open question. 
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APPENDIX A: THE Sp(7V) DIMER MODEL 

As a first step, one introduces bosonic operators 
("Schwinger bosons") {&-[■, &J.}, to represent the spin op- 
erators in the following way: S z = (fejfef — &Jfcj)/2, 

S + = In order to represent SU(2) spins 1/2, one 

needs to supplement the bosonic description with the 
constraint rib = + & = 1 on each site. 

In terms of these operators, the antiferromagnetic 
nearest-neighbour Heisenberg Hamiltonian becomes (J > 
0): 

H = jJ2 S i' S J ( A1 ) 

(v) 

= -./E{( eCTT5 t^/V2)(^^^/V2)-l/4} , 

m 

where Si are spin 1/2 SU(2) operators. This expression is 
formally generalised to Sp(AT) by the introduction of an 
additional flavour index for the bosons, which we label by 
a capital letter. In Sp(JV), there are N such flavours; the 
case of SU(2) corresponds to one flavour: SU(2)~Sp(l). 
The Sp(A^) generalisation, 3j^ Bl of the Levi-Civita sym- 
bol continues to be off-diagonal and antisymmetric in 
the spin index but simply diagonal in the flavour index: 
J ax = €^ v 8ab- We define the singlet (hereafter called 
dimer) annihilation operator for bond (ij): 

en = JZhb£bf T /V2N , (A2) 

where summation over spin and colour indices is implicit. 
The Sp(A^) Hamiltonian reads 

H = - J Y,{ s h^-^}- (A3) 

1. The dimer model at leading order 

A hardcore dimer model is obtained from this Hamil- 
tonian at leading ordeB& as one takes N — > 00 & To see 
this, consider the properties of dimer coverings of the lat- 
tice. Let P denote an ordered pairing of the sites of the 
lattice and \tpp) the corresponding singlet wavefunction: 
\ijjp) = Y\ s l 2n _ lP2n |0), where |0) is the state with no 
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bosons present and P2n-i and p2n are the two members 
of the n th pair contained in P. n runs from 1 to N s /2, 
where N s denotes the number of sites. 

Any such state \ipp) satisfies the constraint rib = 1 for 
every site. At leading order, O(l/N ), these coverings 
are orthogonal, degenerate eigenf unctions of the Sp(A r ) 
Hamiltonian: H\ipp) = Eq \ij>p}, with the ground state 
energy Eq = —JN s /2. Due to the constraint on the 
boson number, it is not possible to have more than one 
dimer per site, and due to the form of the Hamiltonian, it 
is disadvantageous to have less. At O(l/N ), the ground 
states are thus all the hardcore dimer coverings, denoted 
by P, of the lattice. 

2. Derivation of the quantum dimer model 

As the different dimer coverings are not exactly orthog- 
onal, one obtains an overlap matrix, S, between different 
coverings which has the following schematic form: 

S = 1 -en + e 2 o + e 2 n x n + 0(e 3 ) . (A4) 

Here, J is the unit matrix, and the symbols □, o denote a 
nonvanishing matrix element between two dimer configu- 
rations differing in the positions of two and three dimers, 
respectively. 

The signs of the overlap matrix elements cannot be 
chosen freely; in fact, there is no convention which uni- 
formly yields a positive sign in front of the □ term. This 
is in contrast to the case of the square and triangular 
lattices, where this choice is a matter of convention»ii*i^ 



However, one is free to choose the sign of the O terms to 
be, for example, all positive or all negative. 

We can now formally orthonormalise the basis set 
by introducing basis states |\& p ) = J2 P > {^~ 1 ^ 2 ) pp ' IVy)- 
Here, S~ x / 2 is the (symmetric) inverse square root of S, 
which can be obtained from S by a Taylor expansion in s. 
The \^p) can be labelled uniquely by a parent dimer con- 
figuration P, as the other dimer states are admixed only 
to higher order in e. Explicitly, orthonormality follows 
from 

(%\%) = E^X^^vivy) 

p'q' 

= (S-^SS- 1 / 2 )^ = l m . (A5) 

The matrix elements of H in the orthogonalised basis 
now read 

H pq = (%\H\%) =J2s;i /2 S m V 2 (^\H\^ g ,) 
p'q' 

= (S-^HS- 1 ^ . (A6) 

This expression is useful as it is also possible to ex- 
pand (Vy | H \ipq>) in powers of e. It turns out to be con- 
venient to subtract off the energy of a dimerised state, 
— J(l — l/4iV 2 )iV s /2, so that diagonal terms are absent: 
(tp p r | H \ipp>) = 0(e 2 ), and the expansion of the Hamilto- 
nian matrix contains no terms of 0(e ). Thence, 

H QDM = S- 1/2 HS- 1/2 = 2Jen + 0(e 2 ). (A7) 
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